Single nucleus transcriptomics of ventral midbrain identifies glial activation associated with chronic opioid use disorder.

Dynamic interactions of neurons and glia in the ventral midbrain mediate reward and addiction behavior. We studied gene expression in 212,713 ventral midbrain single nuclei from 95 individuals with history of opioid misuse, and individuals without drug exposure. Chronic exposure to opioids was not associated with change in proportions of glial and neuronal subtypes, however glial transcriptomes were broadly altered, involving 9.5 - 6.2% of expressed genes within microglia, oligodendrocytes, and astrocytes. Genes associated with activation of the immune response including interferon, NFkB signaling, and cell motility pathways were upregulated, contrasting with down-regulated expression of synaptic signaling and plasticity genes in ventral midbrain non-dopaminergic neurons. Ventral midbrain transcriptomic reprogramming in the context of chronic opioid exposure included 325 genes that previous genome-wide studies had linked to risk of substance use traits in the broader population, thereby pointing to heritable risk architectures in the genomic organization of the brain's reward circuitry.

Single Nucleus Transcriptomics Reveals Pervasive Glial Activation in Opioid Overdose Cases.

Dynamic interactions of neurons and glia in the ventral midbrain (VM) mediate reward and addiction behavior. We studied gene expression in 212,713 VM single nuclei from 95 human opioid overdose cases and drug-free controls. Chronic exposure to opioids left numerical proportions of VM glial and neuronal subtypes unaltered, while broadly affecting glial transcriptomes, involving 9.5 - 6.2% of expressed genes within microglia, oligodendrocytes, and astrocytes, with prominent activation of the immune response including interferon, NFkB signaling, and cell motility pathways, sharply contrasting with down-regulated expression of synaptic signaling and plasticity genes in VM non-dopaminergic neurons. VM transcriptomic reprogramming in the context of opioid exposure and overdose included 325 genes with genetic variation linked to substance use traits in the broader population, thereby pointing to heritable risk architectures in the genomic organization of the brain's reward circuitry.

Induction of dopaminergic neurons for neuronal subtype-specific modeling of psychiatric disease risk.

Dopaminergic neurons are critical to movement, mood, addiction, and stress. Current techniques for generating dopaminergic neurons from human induced pluripotent stem cells (hiPSCs) yield heterogenous cell populations with variable purity and inconsistent reproducibility between donors, hiPSC clones, and experiments. Here, we report the rapid (5 weeks) and efficient (~90%) induction of induced dopaminergic neurons (iDANs) through transient overexpression of lineage-promoting transcription factors combined with stringent selection across five donors. We observe maturation-dependent increase in dopamine synthesis and electrophysiological properties consistent with midbrain dopaminergic neuron identity, such as slow-rising after- hyperpolarization potentials, an action potential duration of ~3 ms, tonic sub-threshold oscillatory activity, and spontaneous burst firing at a frequency of ~1.0-1.75 Hz. Transcriptome analysis reveals robust expression of genes involved in fetal midbrain dopaminergic neuron identity. Specifically expressed genes in iDANs, as well as those from isogenic induced GABAergic and glutamatergic neurons, were enriched in loci conferring heritability for cannabis use disorder, schizophrenia, and bipolar disorder; however, each neuronal subtype demonstrated subtype-specific heritability enrichments in biologically relevant pathways, and iDANs alone were uniquely enriched in autism spectrum disorder risk loci. Therefore, iDANs provide a critical tool for modeling midbrain dopaminergic neuron development and dysfunction in psychiatric disease.

HIV integration in the human brain is linked to microglial activation and 3D genome remodeling.

To explore genome organization and function in the HIV-infected brain, we applied single-nuclei transcriptomics, cell-type-specific chromosomal conformation mapping, and viral integration site sequencing (IS-seq) to frontal cortex from individuals with encephalitis (HIVE) and without (HIV+). Derepressive changes in 3D genomic compartment structures in HIVE microglia were linked to the transcriptional activation of interferon (IFN) signaling and cell migratory pathways, while transcriptional downregulation and repressive compartmentalization of neuronal health and signaling genes occurred in both HIVE and HIV+ microglia. IS-seq recovered 1,221 brain integration sites showing distinct genomic patterns compared with peripheral lymphocytes, with enrichment for sequences newly mobilized into a permissive chromatin environment after infection. Viral transcription occurred in a subset of highly activated microglia comprising 0.33% of all nuclei in HIVE brain. Our findings point to disrupted microglia-neuronal interactions in HIV and link retroviral integration to remodeling of the microglial 3D genome during infection.

Chromatin architecture in addiction circuitry identifies risk genes and potential biological mechanisms underlying cigarette smoking and alcohol use traits.

Cigarette smoking and alcohol use are among the most prevalent substances used worldwide and account for a substantial proportion of preventable morbidity and mortality, underscoring the public health significance of understanding their etiology. Genome-wide association studies (GWAS) have successfully identified genetic variants associated with cigarette smoking and alcohol use traits. However, the vast majority of risk variants reside in non-coding regions of the genome, and their target genes and neurobiological mechanisms are unknown. Chromosomal conformation mappings can address this knowledge gap by charting the interaction profiles of risk-associated regulatory variants with target genes. To investigate the functional impact of common variants associated with cigarette smoking and alcohol use traits, we applied Hi-C coupled MAGMA (H-MAGMA) built upon cortical and newly generated midbrain dopaminergic neuronal Hi-C datasets to GWAS summary statistics of nicotine dependence, cigarettes per day, problematic alcohol use, and drinks per week. The identified risk genes mapped to key pathways associated with cigarette smoking and alcohol use traits, including drug metabolic processes and neuronal apoptosis. Risk genes were highly expressed in cortical glutamatergic, midbrain dopaminergic, GABAergic, and serotonergic neurons, suggesting them as relevant cell types in understanding the mechanisms by which genetic risk factors influence cigarette smoking and alcohol use. Lastly, we identified pleiotropic genes between cigarette smoking and alcohol use traits under the assumption that they may reveal substance-agnostic, shared neurobiological mechanisms of addiction. The number of pleiotropic genes was ~26-fold higher in dopaminergic neurons than in cortical neurons, emphasizing the critical role of ascending dopaminergic pathways in mediating general addiction phenotypes. Collectively, brain region- and neuronal subtype-specific 3D genome architecture helps refine neurobiological hypotheses for smoking, alcohol, and general addiction phenotypes by linking genetic risk factors to their target genes.

Convergence of case-specific epigenetic alterations identify a confluence of genetic vulnerabilities tied to opioid overdose.

Opioid use disorder is a highly heterogeneous disease driven by a variety of genetic and environmental risk factors which have yet to be fully elucidated. Opioid overdose, the most severe outcome of opioid use disorder, remains the leading cause of accidental death in the United States. We interrogated the effects of opioid overdose on the brain using ChIP-seq to quantify patterns of H3K27 acetylation in dorsolateral prefrontal cortical neurons isolated from 51 opioid-overdose cases and 51 accidental death controls. Among opioid cases, we observed global hypoacetylation and identified 388 putative enhancers consistently depleted for H3K27ac. Machine learning on H3K27ac patterns predicted case-control status with high accuracy. We focused on case-specific regulatory alterations, revealing 81,399 hypoacetylation events, uncovering vast inter-patient heterogeneity. We developed a strategy to decode this heterogeneity based on convergence analysis, which leveraged promoter-capture Hi-C to identify five genes over-burdened by alterations in their regulatory network or "plexus": ASTN2, KCNMA1, DUSP4, GABBR2, ENOX1. These convergent loci are enriched for opioid use disorder risk genes and heritability for generalized anxiety, number of sexual partners, and years of education. Overall, our multi-pronged approach uncovers neurobiological aspects of opioid use disorder and captures genetic and environmental factors perpetuating the opioid epidemic.

Neuron-specific chromosomal megadomain organization is adaptive to recent retrotransposon expansions.

Regulatory mechanisms associated with repeat-rich sequences and chromosomal conformations in mature neurons remain unexplored. Here, we map cell-type specific chromatin domain organization in adult mouse cerebral cortex and report strong enrichment of Endogenous Retrovirus 2 (ERV2) repeat sequences in the neuron-specific heterochromatic B2NeuN+ megabase-scaling subcompartment. Single molecule long-read sequencing and comparative Hi-C chromosomal contact mapping in wild-derived SPRET/EiJ (Mus spretus) and laboratory inbred C57BL/6J (Mus musculus) reveal neuronal reconfigurations tracking recent ERV2 expansions in the murine germline, with significantly higher B2NeuN+ contact frequencies at sites with ongoing insertions in Mus musculus. Neuronal ablation of the retrotransposon silencer Kmt1e/Setdb1 triggers B2NeuN+ disintegration and rewiring with open chromatin domains enriched for cellular stress response genes, along with severe neuroinflammation and proviral assembly with infiltration of dendrites . We conclude that neuronal megabase-scale chromosomal architectures include an evolutionarily adaptive heterochromatic organization which, upon perturbation, results in transcriptional dysregulation and unleashes ERV2 proviruses with strong neuronal tropism.

Common Genetic Variation in Humans Impacts In Vitro Susceptibility to SARS-CoV-2 Infection.

The host response to SARS-CoV-2, the etiologic agent of the COVID-19 pandemic, demonstrates significant interindividual variability. In addition to showing more disease in males, the elderly, and individuals with underlying comorbidities, SARS-CoV-2 can seemingly afflict healthy individuals with profound clinical complications. We hypothesize that, in addition to viral load and host antibody repertoire, host genetic variants influence vulnerability to infection. Here we apply human induced pluripotent stem cell (hiPSC)-based models and CRISPR engineering to explore the host genetics of SARS-CoV-2. We demonstrate that a single-nucleotide polymorphism (rs4702), common in the population and located in the 3' UTR of the protease FURIN, influences alveolar and neuron infection by SARS-CoV-2 in vitro. Thus, we provide a proof-of-principle finding that common genetic variation can have an impact on viral infection and thus contribute to clinical heterogeneity in COVID-19. Ongoing genetic studies will help to identify high-risk individuals, predict clinical complications, and facilitate the discovery of drugs.

Common genetic variation in humans impacts in vitro susceptibility to SARS-CoV-2 infection.

The host response to SARS-CoV-2, the etiologic agent of the COVID-19 pandemic, demonstrates significant inter-individual variability. In addition to showing more disease in males, the elderly, and individuals with underlying comorbidities, SARS-CoV-2 can seemingly render healthy individuals with profound clinical complications. We hypothesize that, in addition to viral load and host antibody repertoire, host genetic variants also impact vulnerability to infection. Here we apply human induced pluripotent stem cell (hiPSC)-based models and CRISPR-engineering to explore the host genetics of SARS-CoV-2. We demonstrate that a single nucleotide polymorphism (rs4702), common in the population at large, and located in the 3'UTR of the protease FURIN, impacts alveolar and neuron infection by SARS-CoV-2 in vitro. Thus, we provide a proof-of-principle finding that common genetic variation can impact viral infection, and thus contribute to clinical heterogeneity in SARS-CoV-2. Ongoing genetic studies will help to better identify high-risk individuals, predict clinical complications, and facilitate the discovery of drugs that might treat disease.

A chromosomal connectome for psychiatric and metabolic risk variants in adult dopaminergic neurons.

BACKGROUND: Midbrain dopaminergic neurons (MDN) represent 0.0005% of the brain's neuronal population and mediate cognition, food intake, and metabolism. MDN are also posited to underlay the neurobiological dysfunction of schizophrenia (SCZ), a severe neuropsychiatric disorder that is characterized by psychosis as well as multifactorial medical co-morbidities, including metabolic disease, contributing to markedly increased morbidity and mortality. Paradoxically, however, the genetic risk sequences of psychosis and traits associated with metabolic disease, such as body mass, show very limited overlap. METHODS: We investigated the genomic interaction of SCZ with medical conditions and traits, including body mass index (BMI), by exploring the MDN's "spatial genome," including chromosomal contact landscapes as a critical layer of cell type-specific epigenomic regulation. Low-input Hi-C protocols were applied to 5-10 × 103 dopaminergic and other cell-specific nuclei collected by fluorescence-activated nuclei sorting from the adult human midbrain. RESULTS: The Hi-C-reconstructed MDN spatial genome revealed 11 "Euclidean hot spots" of clustered chromatin domains harboring risk sequences for SCZ and elevated BMI. Inter- and intra-chromosomal contacts interconnecting SCZ and BMI risk sequences showed massive enrichment for brain-specific expression quantitative trait loci (eQTL), with gene ontologies, regulatory motifs and proteomic interactions related to adipogenesis and lipid regulation, dopaminergic neurogenesis and neuronal connectivity, and reward- and addiction-related pathways. CONCLUSIONS: We uncovered shared nuclear topographies of cognitive and metabolic risk variants. More broadly, our PsychENCODE sponsored Hi-C study offers a novel genomic approach for the study of psychiatric and medical co-morbidities constrained by limited overlap of their respective genetic risk architectures on the linear genome.